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ABSTRACT 

Small amounts of star formation in elliptical galaxies are suggested by several results: surprisingly 
young ages from optical line indices, cooling X-ray gas, and mid-IR dust emission. Such star formation 
has previously been difficult to directly detect, but using UV Hubble Space Telescope (HST) Wide 
Field Camera 3 (WFC3) imaging, we have identified individual young stars and star clusters in four 
nearby ellipticals. This technique is orders of magnitude more sensitive than other methods, allowing 
detections of star formation to 10 -5 M Q yr _1 . Ongoing star formation is detected in all galaxies, 
including three ellipticals that have previously exhibited potential signposts of star forming conditions 
(NGC 4636, NGC 4697, and NGC 4374), as well as the typical red and dead NGC 3379. The current 
star formation in our closest targets, where we are most complete, is between 1 — 8 x 10 -5 Mq yr _1 . 
The star formation history was roughly constant from 0.5 — 1.5 Gyr (at 3 — 5 x 10~ 4 M Q yr _1 ), 
but decreased by a factor of several in the past 0.3 Gyr. Most star clusters have a mass between 
10 2 — 10 4 Mq. The specific star formation rates of ~ 10~ 16 yr _1 (at the present day) or ~ 10 -14 yr" 1 
(when averaging over the past Gyr) imply that a fraction 10 -8 of the stellar mass is younger than 
100 Myr and 10~ 5 is younger than 1 Gyr, quantifying the level of frosting of recent star formation over 
the otherwise passive stellar population. There is no obvious correlation between either the presence 
or spatial distribution of postulated star formation indicators and the star formation we detect. 

Subject headings: Galaxies: elliptical and lenticular, cD — Galaxies: star clusters: general — Galaxies: 
star formation — Ultraviolet: galaxies — Ultraviolet: stars 



1. INTRODUCTION 

As a class, normal elliptical galaxies are not actively 
forming stars. Their star formation rate (SFR) must be 
significantly less than 1 M Q yr _1 , or such activity would 
have been identified long ago in the form of prominent 
H(3 lines or an excess of blue light. However, lower levels 
of star formation are suggested by small changes in opti- 
cal line indices, which can be used to date stellar popula- 
tions of galaxies. At constant metallicity, the Balmer line 
indices are sensitive to age, with inferred ages often sur- 
pr isingly young . For e xample, in the sample of ellipticals 
in iTrager et al.l (|2000D . 40% have mean ages less than 6 
Gyr a nd several are younger than 3 Gyr. iDenicolo et al.l 
(2005) obtain similar results with similar methods. This 
implies that many elliptical galaxies were forming stars 
quite recently (z < 0.5), in conflict with other obser- 
vations that indicate little star formation in elli pticals 
since z » 1 (jDaddi et al -1120011 lLabbe et al.|[200l . One 
explanation for this discrepancy is the "frosting" effect, 
whereby a small rate of ongoing star formation contam- 
inates the Balmer lines, making them relatively strong, 
while contributing little to the mass of the galaxy. The 
degree of the frosting could be understood if the rate of 
ongoing star formation were known, resulting in more 
reliable age determinations. 
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Several mechanisms can lead to small amounts of 
star formation. One avenue involves stellar winds from 
the ensemble of old stars in a galaxy, amounting to 
0.1-1 M Q yr" 1 for ~L* ellipticals. Most of this stel- 
lar mass loss is believed to collide with the ambient 
interstellar medium (ISM) and become heated to X- 
ray emitting temperatures t hat are typically 5 x 10 6 K 
(|Mathews fc Brighcnti 2003), although the efficiency of 
this process is unknown an d some gas may remain cool 
(jParriott fc Bregmanl 120081 ). and occasionally may form 
into new stars. Most of the hot ISM in a galaxy has 
a cooling time much less than a Hubble time, and if 
not driven out as a galactic wind (through active galac- 
tic nuclei or supernovae), the gas will radiatively cool. 
Radiatively cooling gas is found in 30 — 40% of ellipti- 
cals, as apparent from the detection of the O VI line in 
Far U ltraviolet Spectroscop ic Explorer (FUSE) observa- 
tions (Bregma n et al.lfeOOSt h implying a cooling rate of 
~ 0.1 — 0.5 M Q yr _1 . Gas emitting O VI is at the peak 
of the cooling curve, so it will cool to the 10 — 10 4 K 
range. This cooled gas, confined to the central region 
(< 1 kpc), is a natural source of material for star forma- 
tion. In addition to the internal recycling of galactic gas, 
there may be infall of material onto a galaxy from other 
smaller galaxies or from ambient group material. 

Although the gaseous content of elliptical galaxies is 
generally dominated by hot (5 x 10 6 K) X-ray emitting 
material with masses of 10 8 - 10 10 M© (jRoberts et al.1 
1991), some evidence for cool gas (< 10 4 K) exists in 
addition to O VI emission. Ha-emitting material at ~ 
10 4 K is present in the central kpc (10— 20") of most ellip- 
ticals (|Mathewslll9"90HCaon et al.ll2000f ). although it con- 
stitutes a relatively small amount of gas (10 4 — 10 5 Mq). 
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While few elliptical galaxies have detectable amounts 
of H I (< 10 7 - 10 8 Mp) or H 2 ([Roberts et all HMD: 
Idi Sereeo Alighieri et alJl200l. 5-10 % were detected at 
60-100 jum bv IRAS (jBregman et al.lll998l ). showing that 
dust emission is occasionally present. Extinction by dust 
lanes is also seen in ~ 50% of ellipt icals, and is generally 
near the center (|Lauer et al.l 120051 ). Evidently, cool gas 
exists in ellipticals and is most common in the central 
region, but we do not know the fate of this material nor 
whether its presence is stable or varies with time. 

A good strategy for quantifying the amount of re- 
cent star formation in normal ellipticals is to work in 
the ultraviolet (UV), where old, red stars that domi- 
nate these galaxies contribute little compared to hot hor- 
izontal branch stars, post-Asy mptotic Giant B ranch (p- 
AGB) stars, and young stars (jO ' Connelll 1 1 9991 ) . Indeed, 
many near-UV (NUV) observations have been made of 
early-type galaxies in an attempt to quantify amounts of 
recent star formation. For instance, using the Wide-Field 
Planetary Ca mera 2 (WFPC2 ) on th e Hubble Space Tele- 
scope (HSTl. lFerreras fc Silkl (|2000D found that most el- 
lipticals in a z — 0.4 cluster must have young stellar 
mass fractions of at least 0.1%, and in many cases 1 
to 10%, to accou nt for their res t fram e 2000A emission. 
At lower redshift, iKavirai et al.1 (|2007| ) studied a sample 
of 2100 early-type galaxies selected systematically from 
the Sloan Digital Sky Survey and found that 30% of 
them had Galaxy Evolution Explorer (GALEX) NUV- 
optical colors that required a young (< 1 Gyr) compo- 
nent, amounting to 1-3% o f their stellar mass (see also 
the similar earlier study bv lYi et aDl2005f ). 

A handful of UV observations have also been made of 
individual stars, though these observations have mostly 
been limited to the nearest gala xies. Perhaps th e 
most comprehensive set comes from lBrown et al.l (2008), 
where the UV population of stars in the nearby elliptical 
M32 (with a distance, d — 800 kpc) was studied via near 
and far-UV data from the Space Tele scope Imaging Spec- 
trograph (STIS) on HST (see also IBrown et al.Tl2000T ) 
and using stellar e volution models. However, because 
IBrown et a l. (2008) focused on determining the UV con- 
tribution from p-AGB stars and other stars at late evolu- 
tionary stages, and because the field of view of STIS was 
small (25 x 25"), no analysis of young stars in M32 was 
performed. This technique of imaging individual stars 
has also been applied to nearby lenticular galaxies; for 
example, NGC 5 1 02, a nearby SO (d = 4.0 Mpc), where 
iDeharveng et al.l (|1997[ ) used the Faint Object Camera 
(FOC) on HST to search for individual young stars as 
point sources and deduced a li mit for the star forma - 
tion rate of 5 x 10~ 4 M yr _1 . IDeharveng et al.l (|1997f ) 
also observed NGC 3115, another nearby SO (d = 9.7 
Mpc), though they detected no stars and concluded that 
the FOC was not useful for detecting these stars at 
the distances o f most elliptical galaxies. More recently, 
ICrockett et al.1 (|2011l ) detected clumps of blue stars in 
the more distant SO NGC 4150 (d = 13.7 Mpc), using 
HST's Wide Field Camera 3 (WFC3), and determined 
the star formation history of NGC 4150 by fitting the 
UV-optical SED on a pixel-by-pixel basis to a two-burst 
model. 

In this paper we present HST WFC3 UV imaging of 
four nearby elliptical galaxies, which we used to deter- 



mine star formation rates through the direct identifica- 
tion of young stars and star clusters. Star formation 
rates were measured to levels of 10 -5 M Q yr _1 , a limit 
that is much lower than the < 0.01 M Q yr -1 limits that 
have pr eviously been at tained for galaxies at similar dis- 
tances (jO'Connelll 11999). Our sample, which we describe 
in W2.ll includes three ellipticals with potential signposts 
of star forming conditions, plus one control galaxy with 
no prior indication of star formation. Observations and 
photometry are also described in fj2j We present the de- 
tected UV bright sources in Sj3j and discuss their physical 
nature in Sj4l Star formation rates and histories from the 
UV bright clusters are presented in ^ and we conclude 
with a discussion in fj6] 

2. DATA 
2.1. Sample Selection 

Our sample consists of nearby normal ^L* ellipti- 
cal galaxies, which are located between 10.6 Mpc and 
the Virgo Cluster, with the farthest targeted galaxy at 
d = 18.4 Mpc (all distanc es were determine d via surface 
brightness fluctuations bv lTonrv et al.l[2001l ). Three nor- 
mal elliptical galaxies that demonstrated potential sign- 
posts of star formation were included, as well as one typ- 
ical elliptical that has shown no indication of having cur- 
rent star formation. The observed targets are summa- 
rized in Table Q] and are presented below: 

NGC 3379 — We s elected the nearest normal elliptical 
that is old (« 9 Gvr; lTrager et al.l 12000) and has no in- 
dication of recent star formation, NGC 3379 (M105). It 
is the closest galaxy in our sample at d — 10.6 Mpc, and 
was chosen as the control galaxy to allow tests of whether 
the possible star formation indicators in the other galax- 
ies within our sample, such as O VI emission, AGN ac- 
tivity, and PAH emission, are reliable tracers of ongoing 
star formation. 

NGC 4697 — With a distance of 11.8 Mpc, NGC 4697 
is our second closest target and is usually designated a 
typical elliptical galaxy. However, from a sample of 30 
nearby nor mal ellipticals observed with the Spitzer Space 
Telescope (Bregman et al. 2006a, b), it is the only galaxy 
that exhibited a polycyclic aromatic hydrocarbon (PAH) 
feature. Although PAH emission can be produced with- 
out star formation, e.g. , via the winds o f intermediate- 
age stellar populations (V ega et al.l [20l61). it is often an 
indicator of such activity. 

NGC 4636 — A cooling flow elliptical within the Virgo 
Cluster, NGC 4636 has a distance of 14.7 Mpc. It 
has abundant X-ray emitting gas that cools on short 
times cales. and exhibits O VI emission (Bre gman et al.l 
l200ll 120051) . This O VI emission is an excellent tracer of 
cooler gas that could feed star formation and indicates a 
cooling rate of 0.3 M Q yr _1 . 

NGC 4374 — The farthest of our targets, NGC 4374 
(M84) is in the Virgo Cluster and has a distance of 
18.4 Mpc. O VI emiss ion has been detected in this galaxy 
(jBregman et al. 2005), suggesting the presence of cooling 
gas that can feed star formation. Radio jets and AGN 
activity are also pr esent, which may cor relate with star 
formation (e.g., see lHeckman et aT1l2004D . 

2.2. Observations and Data Reduction 

Observations were conducted using HST's WFC3 be- 
tween 2009 November 23 and 2010 March 29 (Cycle 
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Table 1 

Target Summary 



Galaxy 


Type 


m - M 


d 


Rett 


Notes 








(Mpc) 


(arcsec) 




NGC 3379 


El 


30.12 


10.6 


29.9 


no prior star formation indicator 


NGC 4697 


E6 


30.35 


11.8 


42.4 


PAH emission 


NGC 4636 


EO-1 


30.83 


14.7 


59.3 


O VI emission, cooling flows 


NGC 4374 


El 


31.32 


18.4 


34.8 


O VI emission, AGN activity 



N ote. — Galaxy typ es are from RC3 (de Vaucoulcurs ct al. 1991), distances are 
from lTonrv et al.l 1120011). and effective radii are from The 2MASS Large Galaxy Atlas 
(LGA; IJarrett et alj|2003ft . 



17, PID 11583). These observations were made using 
UVIS through two filters, F225W and F336W, which are 
spaced far enough in wavelength that stars on the main 
sequence would stand out in color from possible contam- 
inants such as UV-bright globular clusters on a color- 
magnitude diagram (CMD). Not only are the background 
levels in the UV substantially lower than at optical wave- 
lengths, but fluctua tions in the ba ckground due to hori- 
zontal branch stars ( Worthcy 1993) arc much fainter than 
young stars, which are expected to emit strongly in the 
UV. 

Each galaxy was observed over two orbits for each of 
the F225W and F336W filters. Dithered observations 
were performed over each two-orbit sequence, using the 
UVIS box dithering pattern with four pointings per orbit 
to aid in the removal of artifacts. Each pointing had an 
exposure time of 600s, for a total of 4800s per filter per 
galaxy. As the exposure times were the same for all four 
of our targets, the data of NGC 3379 and NGC 4697, our 
closest targets, are deeper than those of NGC 4636 and 
NGC 4374. 

Data were reduced on-the-fly during retrieval from the 
Multimission Archive at STScI (MAST), resulting in 
pipeline-calibrated and flat-fielded individual exposures. 
We then removed cosmic rays, corrected geometric dis- 
tortions, and combined th e dithered pointings fro m each 
orbit using Multidrizzle (Koekemo er et all 120 02). The 
resulting drizzled images were used as the reference frame 
while the photometry was performed. 

2.3. Photometry 

We per formed PSF p hotometry using DOLPHOT (ver- 
sion 2.0; iDolphinl [2000) . a stellar photometry package 
with specific HST instrument modules. In addition to its 
customized settings for WFC3 data, DOLPHOT can per- 
form localized background subtraction during the pho- 
tometering, which was crucial due to the background 
emission near the center of each galaxy that needed to 
be subtracted carefully, particularly in F336W. 

For each galaxy, we used the deep multidrizzlcd F225W 
image as the reference frame. Initial estimates for the po- 
sition shifts, magnification and rotation were made using 
the wfc3fttdistort task, and the individual dithered expo- 
sures were then photomctered on a chip by chip basis. All 
recommended settings for WFC3 data were used, with 
the exception of the sky-fitting setting (FitSky). Instead 
of fitting the sky prior to each photometry measurement, 
the sky was fit inside the PSF region of 10 pixels but 
outside the photometry aperture of 4 pixels. This was 
done to properly take into the account the varying back- 
ground levels near the center of each galaxy, particularly 
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Figure 1. 80% completeness limits as a function of radius for 
one of our targeted galaxies, NGC 4697, for the F225W (blue) and 
F336W (red) niters. Chips are distinguished by solid (chip 1) and 
dashed (chip 2) lines. While there is a radial trend in the ability 
to detect sources, it is only in the innermost region of the galaxy 
that the variation is significant. 

at radii dominated by background levels in F336W. The 
output magnitudes from dolphot are automatically cal- 
ibrated using Vega zeropoints, while aperture corrections 
were determined manually using the qphot task in IRAF. 
Reddening co r rectio ns are based on the dust maps of 
ISchleeel et afl (fl998l) using R v = 3.fl 

Completeness limits were determined based on artifi- 
cial star tests using approximately 7 x 10 5 artificial stars. 
We imposed a 1/r density distribution for the artificial 
stars, avoiding the chip gap, so that there are approxi- 
mately equal number of stars in uniformly-spaced radial 
bins. Completeness curves were then computed within 
circular annuli centered on each galaxy. We imposed 
magnitude cuts at the 80% completeness level in both of 
the F225W and F336W filters, and then derived the ac- 
curacy of the recovery based on distance from the center 
of the galaxy per chip, as there is a difference in sensitiv- 
ity per UVIS chip. The 80% completeness limit changes 
as a function of radius, with sources farther from the 
center of the galaxy being easier to detect due to the 
lower background. This trend is demonstrated in Fig- 
ure [H which is the 80% completeness limit as a function 
of radius for NGC 4697, per filter and chip. The radial 
trend is weak outside the inner 200 pixels, with only a 
0.3 magnitude maximum variation, where 200 pixels cor- 
responds to only 3.5% of the area within the effective 
radius of NGC 4697. 

After completeness limits were applied, further cuts 
were made to the detected source list based on output 

4 http:/ /wwwmacho.mcmaster.ca/JAVA/Acurve.html 



4 



Ford & Bregman 



from dolphot. First, elongated objects, objects that 
are too sharp, and extended objects were omitted from 
the list (i.e., sources with object types greater than 2). 
All sources containing too many bad or saturated pixels 
were also omitted (i.e., those with error flags greater than 
3). Last, we applied a sharpness cut of |sharpness F22 5w + 
sharpness F336W < 0.4, which we determined via testing 
to be the most appropriate cut for selecting point sources. 

Despite performing the photometry carefully, impos- 
ing strict completeness limits based on the artificial star 
tests, and applying careful cuts based on dolphot qual- 
ity flags, some contamination from extended and back- 
ground sources was apparent in all four galaxy fields as 
well as all eight blank fields that were used to estimate 
background contamination ( £12.41) . We visually identified 
these obvious features, which stood out particularly in 
F336W, and then removed them by hand. We did not 
remove anything that was ambiguous, i.e., we only re- 
moved sources that were clearly elongated and/or clearly 
part of a nearby background galaxy whose bright, star 
forming regions were detected as individual clumps. Af- 
ter removing these features, the remaining sources do not 
appear to correspond to any visible galaxies when their 
positions are overlaid on optical images. 

2.4. Background Control Field 

While many background galaxies are easily identifiable 
by their extended structures, it is certain that unresolved 
background galaxies and point-like AGNs contaminate 
our source lists. To estimate the effect of this contami- 
nation in our fields, we searched MAST for WFC3 obser- 
vations in the F225W and F336W filters that had similar 
exposure times to our data, and used all available data 
that are suitable as blank fields; eight fields from HST 
PID 11359 (Cycle 17; HUDFs 1-8). We retrieved the 
pipeline-calibrated WFC3 F225W and F336W data for 
these fields from MAST and photometered them identi- 
cally to our four galaxies. 

Although the blank fields were observed with the same 
filters and instrument, the exposure times differ from 
those of our observations: each F225W blank field was 
observed for 5688s (roughly 800s longer than our expo- 
sures) and each F336W blank field was observed for 2778s 
(roughly 2000s shorter than our exposures). The similar 
exposure times of the F225W images in the blank fields 
and in our fields, which were used as reference frames in 
all cases, lends to the photometry being performed quite 
similarly with dolphot. As we are particularly inter- 
ested in blue sources, the depth in the redder filter is not 
as important, though it isn't particularly shallow, being 
only 0.3 magnitudes more shallow in F336W. 

Identical cuts based on output from dolphot were ap- 
plied to the blank fields, including object type, error flag, 
and sharpness. Likewise, all obvious extended sources 
that were not removed by these cuts were removed by 
hand, as was done in each of our target fields. We also 
imposed a magnitude cut based on the artificial star tests 
that were run on our galaxies, to ensure that we did 
not include any objects that were only visible due to the 
deeper F225W data. Apertur e an d reddening corrections 
were applied as described in WI.'Sl 

We assume all remaining sources within the blank fields 
to be representative of the unresolved background galaxy 
and point-like AGN population, as there is no reason to 



expect other contaminants such as stars in these fields. 
We combined the source lists of all eight blank fields and 
scaled their total area to the area of one WFC3 field to 
estimate the magnitudes and quantity of expected back- 
ground sources in one WFC3 field. While cosmic vari- 
ance may play a role in the amount of background con- 
tamination, we cannot account for these effects due to 
the small area of the sky covered by the eight fields in 
the HUDF and the lack of similar data in other blank 
fields. 

3. UV-BPJGHT SOURCES 

UV-bright sources are detected in all four of our tar- 
geted galaxies, though there is variation between the 
number and positions of sources per galaxy. An exam- 
ple of several of these sources can be seen in Figure [3J 
which displays a small portion of the observed NGC 3379 
field in both the F225W and F336W filters. Figure [2] 
also demonstrates a noticeable difference in the surface 
brightness of NGC 3379 per filter, where the surface 
brightness in F225W is extremely low, even compared 
to that of F336W. This demonstrates the importance of 
removing the background emission surrounding individ- 
ual sources carefully during the photometry process, as 
discussed in §2.31 It also demonstrates the usefulness 
of the F225W and F336W filters in searching for young 
stars, as the optical filters would be even more domi- 
nated by each galaxy's surface brightness, rendering this 
experiment impossible. 

3.1. Color-Magnitude Diagrams 

Color-magnitude diagrams for all four targeted galax- 
ies are presented in Figure [3l Magnitude limits were 
determined via artificial star tests, and the deepest mag- 
nitude limits are indicated on the CMDs by dashed lines. 
The apparent magnitude limits, which are a function of 
radius, typically ranged from 25.4 < F225W < 25.7 and 
25.5 < F336W < 26.3. These variations were also de- 
pendent on chip, as the UVIS chip 2 is more sensitive 
than chip 1 by 0.2 magnitudes. The inner radial bin, 
with a radius of 100 pixels, had a much shallower lim- 
iting magnitude than the other radial bins, being about 
1 magnitude shallower than the deepest limit, except in 
NGC 4374 where there is a lot of structure in the inner- 
most bin, resulting in the inner bin being about 2 mag- 
nitudes shallower. However, the vast majority of sources 
are not within the first radial bin as it encloses very little 
area (~ 0.2%) of the entire WFC3 field. The absolute 
magnitude limit is strongly dependent on the distance to 
each of these galaxies, resulting in much deeper data for 
NGC 3379 and NGC 4697 than for our other two tar- 
gets, which accounts for much of the discrepancy in the 
number of detected sources for each galaxy. 

Sources are overlaid on a grayscale Hess diagram, 
which indicates the probability of stars at a given loca- 
tion on the CMD assuming a Salpeter initial mass func- 
tion (IMF) and a constant star formation rate (SFR) over 
the past 100 Myr. The Hess diagra m is based on stel- 
lar e volutionary tracks from BaSTI (jPietrinferni et aLI 
l200l . assuming solar metallicity, solar alpha abundance, 
r\ = 0.4 canonical models, and a Salpeter IMF from 2 to 
30 Mq (we are limited at the high mass end by the ab- 
sence of more massive BaSTI tracks, and lower mass stars 
will not be observable with our data). The numbers have 
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Figure 2. A small portion of the smoothed HST WFC3 F225W (left) and F336W (right) images of NGC 3379. Several UV-bright sources 
are detected in both filters (circles). For scale, the white, horizontal line in the bottom, right corner of each image indicates 1" (the entire 
WFC3 field is 162 x 162"). 
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Figure 3. CMDs of targeted galaxies, showing all detected sources with F225W — F336W < 1 (diamonds). Apparent magnitude is on 
the right y-axis while absolute magnitude is on the left y-axis. The grayscale Hess diagram indicates the probability of a star being at a 
given location on the CMD assuming a Salpeter IMF and constant SFR over the pa st 100 Myr, in units of 10 — 6 stars per bin of width 
0.02 mag and height 0.1 mag, and is based on stellar evolutionary tracks from BaSTI (Pictrinfcrni ct al. 2004). The choppy appearance of 
the Hess diagram is an artifact due to the interpolation between different mass tracks. All sources with colors > 1 are omitted from these 
CMDs as they are predominantly globular clusters, which are beyond the scope of this study. Error bars indicate the mean uncertainty, 
as estimated by DOLPHOT, for sources within bins of 0.5 in F225W magnitude. Sources are detected in all of our targets, though there are 
more detected in the nearest targets. Dashed lines show the deepest apparent magnitude limits. 
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b een rescaled to account for lower-mass stars assuming 
a lChabrierl (12001) IMF. 

As globular clusters emit in the UV and are commonly 
found in ellipticals, we determined the colors of globular 
clu sters in the F225W a nd F336W filters using the exist- 
ing [DirscK^E^D (|2005| ) catalogue for NGC 4636, where 
we matched the positions of optically identified globular 
clusters to detected sources in our sample. Out of ~ 80 
matched globular clusters, only 1 had a color blueward 
of F225W - F336W < 1. We therefore expect very lit- 
tle, if any, contamination from globular clusters within 
the color range of F225W - F336W < 1. Redward of 
F225W-F336W = 1, however, the CMDs are domi- 
nated by globular clusters, and so we limit the CMD 
color range to F225W — F336W < 1 as the globular clus- 
ter populations of these galaxies are beyond the scope of 
this study. 

3.2. Source Distribution 

WFC3 F225W images for each galaxy are shown in 
Figure [H where the positions of sources in our sam- 
ple (all sources with F225W - F336W < 1) are marked 
by colored dots. Each color represents a different 
F225W— F336W range: blue represents all sources that 
are near the main sequence (—1.3 < F225W — F336W < 
—0.7), yellow represents all sources with colors to the 
red of the main sequence (F225W - F336W > -0.7), 
and green represents all sources to the blue of the main 
sequence (F225W - F336W < -1.3); there are only 3 
green points, all of which are in NGC 4374. Many sources 
are concentrated around the centers of NGC 3379, 
NGC 4697, and NGC 4374, while there are fewer sources 
near the center of NGC 4636. Also, most sources are not 
within the color range of the main sequence, and instead 
have colors within -0.7 < F225W - F336W < 1. 

3.3. Radial Surface Densities 

Radial surface density profiles of the detected sources 
are given as a function of effective radius in Figure O 
The dotted, horizontal lines in each plot represent the 
mean density of background sources determined from 
the background control field, the CMD for which is 
shown in Figure [5] and was determined by combining the 
eight blank HUDF fields. This background has not been 
subtracted from the target galaxies CMDs in Figure [3l 
The solid, curved lines represent de Vaucouleurs profiles 
(|de Vaucouleurslll948l ) that are scaled by each galaxy's 
corresponding i? c ff and are arbitrarily scaled in ampli- 
tude for reference. Different F225W— F336W ranges are 
represented by different colors, where blue represents 
sources that have colors expected of individual young 
stars based on stellar evolutionary tracks, green spans 
colors expected of background galaxies and AGNs (see 
Figure O, as well as other potential objects that are 
off the main sequence, and red represents the remaining 
sources with F225W— F336W redder than that expected 
from the background sample and main sequence stars. 
The surface densities account for unobservable regions 
within each radial bin, such as the corners of the fields 
and the chip gap. 

The majority of sources are in the intermediate ("back- 
ground") color range. In both NGC 3379 and NGC 4697, 
the surface density of sources in all three color ranges are 



above the mean density of background sources and reveal 
that these sources are concentrated about the centers of 
these galaxies. These trends are also seen in NGC 4374, 
but at lower significance. There is no clear evidence 
for an excess above the background in NGC 4636, ex- 
cept possibly for the reddest sources at intermediate ra- 
dius. This is in agreement with the impression given 
from Figure [H where the sources do appear concentrated 
towards the center of the galaxy, especially those within 
the intermediate color range. Some fraction of these UV- 
bright sources must therefore be associated with the tar- 
get galaxies and cannot simply be an overlap in back- 
ground sources. 

4. SOURCE IDENTIFICATION 

A variety of sources emit in the UV, including old, 
horizontal branch and p-AGB stars, young stars, he- 
lium burning stars on the blue loop, open and globular 
clusters, AGNs, and background galaxies. It is there- 
fore crucial to determine the expected colors and mag- 
nitudes of each of these sources and address the likeli- 
hood of their detection in our targeted galaxies, which 
we describe in the following subsections. As discussed in 
£13.11 we have already eliminated globular clusters from 
our sample by restricting our analysis to sources with 
F225W-F336W < 1, as sources with F225W-F336W 
greater than 1 are much too red to be the young stars 
and open clusters that we are searching for. Globular 
clusters are therefore omitted from the discussion below. 

4.1. Main Sequence Stars 

We used stellar evo lutionary tracks from BaSTI 
(|Pietrinferni et ail [2004). which were derived specifically 
for WFC3 filters, to construct a Hess diagram assuming a 
Salpeter IMF and constant SFR over the past 100 Myr. 
The tip of the main sequence is clearly visible in Fig- 
ureO spanning -1.1 < F225W - F336W < -0.7, and is 
where young, hot stars would be expected on the CMD if 
they were present in our targets. In all but one of the tar- 
geted galaxies there is at least one detected source within 
this color range, making it entirely possible that they are 
young O or B stars. For reference, main sequence O stars 
range from Mp225W = —8.8 to —6.8, while main sequence 
B stars rapidly become too faint to be detected (e.g., a 
B2V star has M F22 5w = -4.8). 

It is clear in the CMDs, however, that the majority of 
sources are not located on the main sequence, but instead 
are mostly redder in color. To determine whether main 
sequence stars could be found at these colors, we esti- 
mated the effect of metallicity on the location of the main 
sequence. We created Hess diagrams using a-enhanced 
models with Z = 0.0001, 0.01, 0.03, and 0.04. Higher 
metallicities indeed shift the evolutionary tracks to the 
red, but the magnitude of the effect is small: e.g., using 
Z = 0.04 (« 2.5Zq) instead of solar metallicity results 
in a shift of the main sequence of +0.1 in color. 

Another possibility that could account for the appear- 
ance of the sources being redder than that of the expected 
location of the main sequence is the presence of dust. It 
is possible that there is dust at the center of some or all 
of our targeted galaxies, and indeed is evident by obvious 
dust lanes that are visible in our data of NGC 4697 and 
NGC 4374. However, these visible dust lanes are limited 
to the very central regions of the galaxy (within a radius 
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Figure 4. Smoothed HST WFC3 F225W images of observed galaxies. Colored dots represent the positions of UV-bright sources in 
our sample, and are color-coded by their F225W — F336W value: blue dots represent sources that are near the main sequence (—1.3 < 
F225W — F336W < —0.7), yellow dots are sources with colors to the red of the main sequence, and green dots are sources to the blue of 
the main sequence. The dot-dashed red annuli show the effective radius. UV-bright sources are detected in all four galaxies, though there 
is variation between the number and position of sources per galaxy. For example, many sources are concentrated around the centers of 
NGC 3379, NGC 4697, and NGC 4374, while there are very few sources near the center of NGC 4636. The total WFC3 field is 162 X 162". 
For reference, the white, horizontal line in the bottom, right corner denotes 30 ". The sources that are visible on these images without 
overlaid dots did not meet our selection criteria. 



of 7" for NGC 4374, the galaxy with the most prominent 
dust features), and most of our targets are much farther 
from this region. So, although some sources may appear 
redder due to dust, it is unlikely that this is the case for 
all of the targets between -0.7 < F225W-F336W < 1.0. 
We have also correc ted for foreground re ddening based 
on the dust maps of ISchlegel et al.l ([1998D . 

We conclude that if the sources lie on the grayscale 
of the Hess diagram in Figure [3l then by definition the 
sources can be individual stars with masses of < 30 M Q , 
but that none or very few of the redder sources are indi- 
vidual stars. It is possible, however, that they represent a 
population of open star clusters, a theory we investigate 
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4.2. Star Clusters 

If stars are forming in these nearby elliptical galaxies, 
it is likely that 70 — 90% of the s tars form embedde d 
in clusters instead of individually (|Lada fc La da 2003). 
Although zero-age clusters are dominated by the most 
massive star, and therefore would lie on the grayscale in 
Figured! they quickly evolve to redder colors and become 
fainter as they lose their bright O and B stars. 

A CMD for each galaxy is shown in Figure [7] 
with evolutionary tracks of open clusters overlaid. 
The evolutionary tracks were determined assuming a 
Salpeter IMF and using BaSTI stellar evolutionary 
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Figure 5. Radial surface density profiles of UV-bright sources within target ed galaxies (points). H orizontal dotted lines are the mean 
density of background sources and solid curves are de Vaucouleurs profiles (dc Vaucoulcurs 1948) arbitrarily scaled in amplitude for 
reference. Different colors represent different F225W— F336W ranges, where blue spans colors expected to be dominated by main sequence 
stars (—1.3 < F225W — F336W < —0.8), green spans the range expected to be dominated by open clusters and background galaxies 
(-0.8 < F225W - F336W < 0.3), and red spans colors greater than these two ranges (0.3 < F225W - F336W < 1.0). Both the red and 
blue horizontal dotted lines are at very low levels, while the green dotted lines are elevated, showing the background galaxy level. 
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Figure 6. CMD of background control field, determined by com- 
bining eight blank HUDF fields. All sources are confined to colors 
between —1.1 and 0.8, with most between —0.8 and 0.3, suggesting 
that the majority of contamination will be found within this range. 

tracks (|Pietrinferni et aLl l2004f l . as a funtion of cluster 
mass, assuming that stars die instantly once they evolve 
off of the main sequence. The main sequence lifetimes, 
tms, are assumed to follow: 

whi ch is a good approxim ation to the lifetimes in Figure 
4 of lSchaller etaU fl992T) over the 30 Myr < r MS < 2 Gyr 
range of relevance. Tracks for clusters with masses rang- 



ing from 10 2 M to 3 x 10 4 M Q are overlaid. These CMDs 
have also been corrected for background contamination, 
where points on the CMD were matched to the distribu- 
tion of sources in the background CMD, which was scaled 
to the area of one WFC3 field, have been removed for a 
clearer picture of the intrinsic CMD. As this procedure 
is statistical in nature, it should not be used when exam- 
ining individual sources, but it gives an accurate picture 
of the effect of background contamination on the CMD. 
Cluster tracks are the same for all four galaxies. 

If the sources are open clusters, their distribution in the 
CMD reflects the cluster mass function, the rate at which 
clusters traverse the CMD, and the recent star formation 
history (SFH). For a constant SFH, the first two factors 
approximately cancel out at a given absolute magnitude: 
high mass clusters are produced less often than low mass 
clusters, but evolve more slowly (and therefore "pile up" 
on the CMD) at a given luminosity because they are later 
in their evolut ion. For exa mple, if the cluster IMF slope 
is —2 (e.g., see Giclcs 200§), clusters with masses between 

10 2 M Q to 3 x 10 2 M Q are 10 times more numerous as 

10 3 M Q to 3 x 10 3 M clusters, but evolve 14 times 
faster between — 6 < F225W < —5, where most of the 
sources are located. We therefore expect star clusters 
to be approximately uniformly distributed in color at a 
given magnitude, as observed. 

We estimated the disruption timescale of the clusters 
using equation (9) in lLamers et al.l (|2005f ) . which scales 
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Figure 7. CMDs with ev olutionary tracks of clus ters overlaid. Cluster tracks assume a Salpeter IMF and are derived from the BaSTI 
stellar evolutionary tracks HPietrinferni et al.|[200D . Tick marks indicate ages of 0, 100, 200, 400, 600, 800, 1000, 1200, 1400, 1600, 1800, 
and 2000 Myr, while dashed lines indicate 100 Myr (left) and 1 Gyr (right). Magnitude limits restrict the age at which we can detect lower 
mass clusters. These CMDs have been corrected for background contamination. 



as cluster mass to the 0.62 power and as the inverse 
square root of the ambient density. This latter quantity, 
which is dominated by the baryonic component at these 
small radii, we estimate using the Uaffq (|1983[ ) model, 
which provides a good model for elliptical galaxy light 
profiles. For a mass of 2 x 10 10 M Q and an effective ra- 
dius of 2 kpc, typical of the target galaxies, the ambient 
density at the effective radius is ~ 0.05 M pc~ 3 , and the 
cluster disruption timescale is > 1 Gyr for 10 4 M Q clus- 
ters, and > 100 Myr for even 10 2 M clusters. We would 
therefore expect that clusters from present-day star for- 
mation have not yet disrupted. Moreover, disruption is 
not instantaneous, and clusters in the process of dissolv- 
ing may well still appear as point sources at these dis- 
tances. 

To determine the radial extent of the sources, it was 
necessary to stack all sources to increase the signal-to- 
noise. A variety of stacks were created spanning differ- 
ent color ranges to ensure that there was no dependence 
of the radius on color; color had no effect on the deter- 
mined radii. Stacks were scaled by the inverse of the 
flux, making all sources weighted equally regardless of 
magn itude. We used T iny Tim HST PSF modeling soft- 
ware (|Krist et al.ll2011h to determine the PSF for WFC3 
F225W filter, resulting in a PSF with FWHM = 1.7. 
The stacked radius was more extended than the PSF, 
with a FWHM — 2.5. The radius of the stacked sources 
for the F336W filter was also more extended than the 
PSF derived using Tiny Tim. While the minimum re- 
solvable radius is 3.4 pc, we determine that the mean 



cluster radius is 5.0 pc (based on the radius of the F225W 
cluster stack). This is consistent with the sources be- 
ing open clusters; for example, the median cluster ra- 
dius in the Catalogue of Open Cluster Data (COCD; 
iKharchenko et alJ l2005) is 4.3 pc. 

4.3. Post-Main Sequence Stars 

Post-main sequence stars, such as hot horizontal 
branch and p-AGB stars, contribute much of the UV 
light from galaxies (jO'Connelll [l999l) . Typical p-AGB 
stars detected in M32 and modeled f rom stellar evolu- 
tionary tracks by iBrown et al.l (|2008l ) are too faint by 
several magnitudes relative to our data to contaminate 
our sample. Theoretical considerations also make it un- 
likely that the sources we detect a re old p-AGB stars. 
The p -AGB evolutionary tracks of iVassiliadis &; Wood! 
(1994) have been overlaid on our deepest CMD in Fig- 
ure where the bolometric correction has been com- 
puted assuming a blackbody spectrum. Although the 
colors of our sources are consistent with being p-AGB 
stars, the magnitudes are not: only stars with initial 
masses > 1.5 M Q could be bright enough to be the 
sources we detect, which implies that they would have 
to belong to a reasonably young (< 4 Gyr) population. 
Moreover, such massive stars live very short lifetimes as 
p-AGB stars, as seen by the length between 100 yr tick 
marks on Figure [5) As these stars spend more time blue- 
ward of the main sequence than redward of it, the com- 
plete absence of sources blueward of the main sequence 
is inexplicable if they are p-AGB stars, but is a natural 
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Figure 8. CMD of UV-bright sources in NGC 3379, our closest 
target and therefore deepest CMD, Overlaid curv e s repr esent p- 
AGB evolutionary tracks from Vassiliadis & Wood (1994). Tracks 
are labelled according to the stellar model initial mass, Mi, and 
the final core mass, M c . Ticks are plotted every 100 yr. A main 
sequence turnoff mass of Mj = 1.0 corresponds to a 10 Gyr popula- 
tion, Mi = 1.5 corresponds to an age of < 4 Gyr and brighter tracks 
are produced by turnoff stars from even younger populations. The 
bottom track is therefore the only place we would expect to see a 
population of p-AGB stars from the dominant old population, and 
our sources are too bright to be these p-AGB stars. This CMD has 
been corrected for background contamination. 

consequence if they are star clusters that evolve redward 
from locations slightly to the red of the main sequence. 
Horizontal branch stars are similar in temperature to p- 
AGB stars, but significantly fainter, and so cannot be 
the sources we detect. 

Helium burning stars on the blue loop are expected to 
have colors redder than the main sequence, though that 
number should be low based on the grayscale probability 
distribution derived from stellar evolutionary tracks, i.e., 
the blue loops are the regions seen in Figure [3] between 
-0.5 < F225W - F336W < 0.2 at the magnitudes we are 
detecting. The grayscale probability distribution clearly 
illustrates that few stars are expected in this region, i.e., 
only the helium burners on the blue loop would be in 
these color ranges, and even then they are so short-lived 
that they would be very unlikely. Also, if these were blue 
loop stars, many more stars on the main sequence would 
be expected based on the Hess diagram, i.e., at magni- 
tudes brighter than Mp225W = —4.6, we would expect 
to see 19 times more stars on the main sequence than on 
the blue loop, though we only detect a few sources on the 
main sequence while detecting many sources redward of 
this. 

4.4. Point-like Background Galaxies and AGNs 

Contamination from point-like background galaxies 
and AGNs is unavoidable. However, the amount of con- 
tamination can easily be estimated by using a blank con- 
trol field to get an estimate of the colors, magnitudes, and 
number of sources expected in such a field. In section 
I2.4l we described eight blank fields observed with WFC3 
using the same filters as our observations, which we pho- 
tometered to derive a background sample. After scaling 
the total area of the combined blank fields to the area of 
one WFC3 field, we determined how many background 
sources would be expected on our CMDs and in what po- 
sition on the CMD they would be in. We then matched 
these background sources to the most likely counterpart 
in our CMDs based on their magnitudes and colors, and 



then subtracted them to determine what the CMD cor- 
rected for background contamination would be. The re- 
sulting CMDs of this process can be see in Figure [3 
which demonstrates that although a background sample 
has been removed from the CMD, there are still numer- 
ous sources within the CMD that are unaccounted for, 
in each of our galaxies. Although galaxy clustering can 
change the magnitude of the background counts, the ra- 
dial dependence of the source distribution clearly demon- 
strates that the majority of detected sources are con- 
tained in the target galaxy and are therefore not point- 
like background galaxies or AGNs. An estimate of the 
variance in the background can be obtained by looking 
at the field-to-field dispersion in the eight blank fields; 
the mean number of sources is 11.6, with a standard de- 
viation of 5.8. In contrast, the number of sources de- 
tected within the color range of our study in NGC 3379, 
NGC 4636, NGC 4697, and NGC 4374 are 98, 31, 87, 
and 41, respectively. 

4.5. Surface Brightness Fluctuations 

It is important to address the likelihood of false de- 
tections due to surface brightness fluctuations (SBFs), 
which will not only result in an increased number of 
detections, but can also mimic a radial ly concentrated 
population. Models by iWorthevI ()1993l ) determine the 
characteristic SBF magnitude in F284W and the opti- 
cal U band for a variety of old stellar populations. In 
the case where there is both an extended/blue horizon- 
tal branch and a low-mass p-AGB star population con- 
tributing to the SBFs, the characteristic SBF magni- 
tude will be the luminosity-weighted mean luminosity 
of the populations. Using their low-mass p-AGB case 
as a worst case scenario, (F284W, U) = (1.2, 2.7), which, 
when converted to the filters we are interested in, implies 
(F225W, F336W) = (-0.3,2.7). As the deepest limit in 
our data is -4.5, this corresponds to nearly a 50er fluctu- 
ation, which is extremely unlikely, particularly over the 
WFC3 field of view. We therefore rule out the possibility 
that these sources are surface brightness fluctuations. 

5. RESULTS 

As the detected sources are most likely individual main 
sequence stars and star clusters, we can deduce their 
individual ages and masses from evolutionary tracks. 
The properties of our sample of UV-bright sources for 
each galaxy are presented in Table [3J where the galaxy 
ID is in column 1, right ascension (a) and declina- 
tion (S) are in columns 2 and 3, apparent F225W and 
F336W magnitudes are in columns 4 and 5, and color 
(F225W-F336W) is in column 6. The derived cluster 
ages and masses are also presented in Table [5] in columns 
7 and 8. Errors for the cluster properties were deter- 
mined by Monte Carlo sampling the error ellipse in the 
CMD around each cluster, and represent the 16th and 
84th percentiles of the Monte Carlo distribution (ler). 
Because the clusters evolve faster at younger ages, the 
errors are asymmetric and there is expected to be a net 
bias in the age and mass determinations. We estimate 
the effect of this by examining the mean of the Monte 
Carlo error distribution. We find that the mean age after 
accounting for the photometric error is ~ 30 Myr older 
than the input age for clusters younger than 1 Gyr (the 
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method is not well denned for older clusters because the 
Monte Carlo points run off of the calculated evolutionary 
tracks), and the recovered mean mass is larger than the 
input mass by a factor falling from 3 at 100 M to ~ 1.5 
at 1000 M and is negligible at 10 4 M . This system- 
atic error should be kept in mind when interpreting the 
cluster properties. 

Assuming that all star formation results in an extant 
cluster, we can add together the mass of clusters within 
a given age range to estimate the star formation rate at 
that time. Star formation rate as a function of cluster 
age is shown in Figure [9l where the individual points 
represent the star formation rate contributed by clusters 
as a function of mass and the solid line represents the 
total star formation rate per bin. The cluster age er- 
rors are typically < 100 Myr for ages < 200 Myr, rising 
to 200 — 400 Myr for older clusters. We have therefore 
adjusted the bin sizes in Figure [9] to be similar to the 
typical errors. There is ongoing star formation in all 
four of the observed galaxies, and there is variation be- 
tween galaxies. As the data get shallower, there is a bias 
that lower mass clusters are lost first. We have applied 
a completeness correction, which assumes that the clus- 
ter catalog is complete to the mass where the cluster age 
tracks cross the typical magnitude limit for each galaxy, 
using a cluster initial mass function that is a power law 
of slope —2 from 10 2 to 10 6 M Q 0. As can be seen in Fig- 
ure[9j applying the completeness correction does not have 
a dramatic effect on the implied star formation rates, and 
we are therefore directly constraining the star formation 
rate with the observed clusters. 

The derived SFRs take background contamination into 
account by subtracting those sources that best match 
the distribution of the scaled control background field. 
At ages > 10 9 yr, the colors of star clusters become 
more sensitive to cluster metallicity than to cluster age 
- changes of 0.1 dex in metallicity become enough to 
move the clusters into different age bins — and therefore 
the ages of clusters we list as older than this should be 
considered with care. Similarly, the details of the star 
formation histories beyond 1 Gyr are very uncertain. As 
noted above, photometric errors cause a mean overesti- 
mate in cluster masses at the low mass end, resulting 
in a tendency to overestimate the star formation rate 
contribution by those clusters. However, only for the 
lowest-age bin in NGC 4636 do the lowest-mass clusters 
contribute a significant amount to the total estimated 
star formation, so this should have a minimal effect on 
the derived star formation rates. The mean age overes- 
timate due to photometric errors is significantly smaller 
than the size of the age bins, so it should not have a 
significant effect. 

In three of the four galaxies, we detect current star 
formation at rates that range from 1.0 x 10~ 6 M Q yr -1 
for NGC 4636 to 8.0 x 10" 5 M Q yr" 1 for NGC 3379 
(note that because the completeness is different be- 
tween galaxies, we quote completeness-corrected val- 
ues here, although the difference is minimal for the 
present-day SFR); in NGC 4374, the lack of detected < 
100 Myr clusters corresponds to an upper limit of 1.15 x 

5 The total mass has a very mild logarithmic dependence on the 
upper and lower limits of the mass function, so the choice of these 
particular values does not introduce a large uncertainty. 



10 6 M yr 1 . Star formation within the past Gyr is de- 
tected in all four galaxies, with average rates somewhat 
larger, ranging from 2.9 x 10~ 4 M Q yr" 1 for NGC 4374 
up to 4.6 x 10~ 4 M Q yr" 1 for NGC 4697. Encouragingly, 
these numbers are of the s ame order as the average star 
formation rates found by iDonovan Mever et al.l (2010) 
for early-type galaxies that had no other evidence for 
star formation, based on total GALEX NUV emission. 
It therefore seems likely that the emission they see is the 
unresolved sum of the sources we detect with HST. 

The cluster mass function presents another projection 
of the cluster data and is shown in FigureQIilfor the near- 
est galaxy in our sample, NGC 3379. The data for the 
remaining galaxies are too shallow to draw strong conclu- 
sions on the mass function, but they are consistent with 
what is seen in NGC 3379. Because the mass function 
is complete to different masses at different ages, we have 
plotted the mass function for clusters in three different 
age ranges, with the completeness limit shown in each 
case as the vertical dotted line. The black diagonal line 
shows a power law of slope —2 for reference. 

One notable feature in Figures [9] and [10] is the lack 
of very massive (> 3000 M©), very young (< 300 Myr) 
clusters. Such clusters would not be too faint to be de- 
tected, and they should not have been disrupted. We can 
imagine several potential explanations for this, however 
most are not viable upon further examination. For in- 
stance, particularly luminous sources could have diffrac- 
tion spikes that are artificially broken up by the source 
finding algorithm; however, the successful detection of 
more luminous sources at redder colors refutes this expla- 
nation. Moreover, visual inspection of those few detec- 
tions that dolphot removed for being extended revealed 
all to be truly extended objects rather than diffraction 
spikes. The systematic age overestimate from photomet- 
ric errors also cannot be the explanation, as this would 
be more important for low mass clusters than high mass 
clusters. Another potential explanation is that there are 
no high mass clusters at any age and the apparent old 
high mass clusters are actually blue loop stars. However, 
as argued in Section 14. 3[ the lack of observed main se- 
quence stars precludes a large population of blue loop 
contaminants. Furthermore, such a population would 
produce a peak in Figure [TO] in the 3000 - 10000 M , 
200 — 600 Myr bin, where no such feature is seen. While 
it is true that the cluster evolution tracks are based on 
stellar models that are poorly constrained in the UV, 
and could potentially be systematically too blue result- 
ing in systematically high age estimates, this also cannot 
be the explanation because it would require the highest 
mass stars to be significantly bluer than the lower mass 
stars. 

The most likely explanation is that there is a physical 
reason why galaxies in this sample have star formation 
that is not vigorous enough to generate clusters more 
massive than 3000 M Q . Given that the galaxies are all 
very red ellipticals, this is entirely possible, particularly 
as they have current star formation rates that are ex- 
tremely low even relative to their average over the past 
Gyr. This is augmented by small number statistics - 
at these low star formation rates, only one very massive 
very young cluster per galaxy would be expected. 

Figure [TT] shows how the current and recent star for- 
mation of the galaxies vary as a function of their stel- 
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Table 2 

Properties of UV-Bright Clusters 



Galaxy 


a 


<5 


F225W 


F336W 


F225W-F336W 


Mass 


Age 
















(y) 


NGC 3379 


10:47:43.078 


+12:35:20.59 


25.04 + 0.19 


24.94 ± 0.08 


0.10 + 0.21 


7.041!;?! x 10 3 


7.06j£M x 10 8 




10:47:44.170 


+ 12:34:42.50 


25.37 + 0.25 


26.03 + 0.19 


-0.67 + 0.32 



10:47:44.325 +12:34:34.47 25.03 + 0.19 25.16 + 0.09 -0.13 + 0.21 3.07lf % X 10 3 4.00l^; 8 9 X 10 ; 

10:47:44.539 +12:34:48.99 24.98 + 0.18 25.37 + 0.11 -0.38 + 0.21 1.171J- $j X 10 3 2.191q X 10 1 

10:47:44.621 +12:34:21.47 25.48 + 0.30 26.15 + 0.21 -0.67 + 0.37 

10:47:45.475 +12:35:01.37 24.60 + 0.14 24.44 + 0.05 0.17 + 0.15 1.301°-™ X 10 4 8.23!^ X 10 ; 



Note. — Tabic [2] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. Sources that lie outside the cluster evolution tracks do not have mass or age listed. Sources 
where the upper age error lies outside the cluster evolution tracks arc listed as lower limits. 
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Figure 9. Star formation rate as a function of cluster age. Individual points represent the star formation rate contributed by clusters 
as a function of mass, while the solid line represents the total star formation rate per bin. There is ongoing star formation in all four of 
the observed galaxies, though more in some than in others. These have been corrected for background contamination. The completeness 
correction assumes that the cluster catalog is complete to the mass where the cluster age tracks cross the typical magnitude limit for each 
galaxy, using a cluster initial mass function that is a power law of slope —2 from 10 2 to 10 6 Mq. Error bars are plotted for the "Total" 
line, and are based on bootstrap resampling of a Poisson-distributed number of clusters within each age bin. The relative errors on the 
completeness-corrected SFRs are identical. 



lar mass, wher e the stellar masse s were calculated from 
2MASS LGA (IJarrett et al.ll2003l) K tnt magn itudes and 
the M/L calibration of iBell fc de Jonel ([2001) using the 
V — K color from the RC3 Vt magnitudes. The Gyr- 
averaged SFRs are more or less independent of stellar 
mass, while the current star formation appears to decline 
with mass. Lines of constant specific star formation rate 
(SSFR) are overplotted. The SSFR values range from 
< 8 x 10~ 18 yr -1 to 2 x 10~ 15 yr" 1 at the present day, 
and 2 x 10 -15 yr" 1 to 10~ 14 yr" 1 averaged over the past 
Gyr. 



6. DISCUSSION 

This work presents the first derivation of star forma- 
tion rates of true "red and dead" ellipticals via direct 
detection of point sources in the UV using the WFC3 
on HST. Our targets have distances ranging from 10.6 
Mpc to the Virgo Cluster at 18.4 Mpc, and three were 
chosen because they demonstrated physical evidence of 
possibly harboring young stars, in the form of PAH emis- 
sion, OVI emission, and AGN activity, while the fourth 
was selected as the control galaxy with no star forma- 
tion indicator. The detected sources are best explained 
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Figure 10. Cluster mass function for NGC 3379 for three different 
age ranges, with the completeness limit shown in each case as the 
vertical dotted line. The black diagonal line shows a power law of 
slope —2 for reference. 
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Figure 11. Star formation rate as a function of stellar mass for 
all four galaxies in our sample. Filled symbols indicate the SFR 
over the past 100 Myr, while open symbols are the average SFR 
over the past Gyr. Gray lines of constant specific star formation 
rate are overlaid. The Gyr-aged SFRs are more or less independent 
of stellar mass, while the current star formation appears to decline 
with mass. 

as being individual star clusters with masses between 10 2 
and 10 4 M & and ages less than 1 Gyr: they are radi- 
ally concentrated about the center of each galaxy, rul- 
ing out background contamination; they are too bright 
to be p-AGB stars from the dominant old population 
and too red as a population to be p-AGB stars from an 
intermediate-age population; and the absence of a large 
number of bright main sequence turnoff stars precludes a 
large enough population of blue loop stars to account for 
any but a tiny fraction of the observed sources. Our de- 
tections of young stars and star clusters in all four of our 
targeted ellipticals confirms that low-level star formation 
is ongoing in these "red and dead" galaxies. 

All of the observed galaxies have line index ages 
(iKuntschner et all [20101: iSanchez-Blazquez et all 120061: 
Trag er et all 12000) that are uniformly very old. To es- 
timate whether the very small amount of star formation 
we detect is consistent with these ages, we compared the 
total F336W light coming from our detected sources to 
the total U-band flux from RC3. For NGC 3379, where 
we go deepest down the young cluster mass function, 
only 0.03% of the U-band light within the effective ra- 
dius is coming from our sources. It is therefore con- 
sistent with the old ages determined spectroscopically. 



A current SFR that is 1-2 orders of magnitude larger, 
which would still be immeasurable by other techniques 
but is easily ruled out for the galaxies in this sample, 
would be required to provide enough frosting to begin to 
contribute enough light at the age-sensitive wavelengths 
to alter the derive d ages. Our data a re also consistent 
with the models of iRogers et al.l (|2010[) , which were per- 
formed to study the star formation history of NGC 4697 
(amongst other ellipticals). Their results show that a 
simple stellar population model has an equally good fit 
to the data as their best fit frosting model, meaning that 
up to 3% of the stellar mass in a younger component can 
be added without changing the fit. Our derived value is 
consistent with this result, as the young population we 
detect makes up a much smaller percentage than this. 

The typical star formation rates of these galaxies are ~ 
10~ 5 M yr" 1 . The two closest galaxies, NGC 3379 and 
NGC 4697, where we have the best measurements of the 
SFR, both show current SFR at a similar level, despite 
the presence of PAH emission from one and not the other. 
The errors are much larger for the more distant galaxies 
with OVI emission, NGC 4636 and NGC 4374, but they 
also show consistent results. This perhaps indicates that 
such emission is a poor tracer of the small-scale minute 
amounts of cold gas that must be present to form stars. 

Comparison to the current stellar mass of the galax- 
ies gives specific star formation rates of ~ 10~ 16 yr _1 at 
the present day, which implies that stars younger than 
100 Myr provide a frosting that consists of 10 -8 of the 
total stellar mass. The Gyr-averaged specific star forma- 
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tion rate is somewhat higher, 
ing to 10~ 5 of the stellar mass. 

There is no obvious structure to the spatial distribution 
of clusters beyond their general concentration towards 
the center of each galaxy, or any obvious correlation be- 
tween the locations of the clusters and any other feature 
in the galaxy. 

The star formation histories of the galaxies reveal a 
drop in the star formation within the past 300 Myr, par- 
ticularly in NGC 3379 and NGC 4697 where the data are 
the deepest and the sources are most obviously physically 
associated with the galaxy. There are several possible 
origins for this, such as a contaminant that we have not 
properly accounted for that pollutes the region of the di- 
agram where older clusters lie, though this explanation 
is highly unlikely. Another possibility is that the clus- 
ter evolution tracks evolve too slowly in this region of 
the CMD. More sophisticated population synthesis mod- 
els in the WFC3 UVIS bands would be required to test 
this hypothesis, which are not presently available. Or 
perhaps our selection of galaxies, which was not random 
but specifically targeted old red galaxies with no known 
star formation, picked out galaxies that are unusual qui- 
escent at the present day. The youngest clusters are the 
easiest to see, easier than old clusters, but they are com- 
pletely absent in 2 of the 4 galaxies. One of these galax- 
ies, NGC 4374, hosts an AGN, which may have recently 
heated the gas and made it unavailable for star forma- 
tion. This hypothesis can be tested with a larger and 
more representative sample of elliptical galaxies, which 
we hope to obtain. 

These results provide a key observable for the amount 
of residual star formation in quenched early-type galax- 
ies. By comparing the rate at which stellar mass transi- 
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tions across the "green valley" to the stell ar mass func- 
tion o f the red and blue galaxy populations. iMartin et al.l 
(2007) inferred a quenching timescale on the order of 
hundreds of Myr. However, because these small levels 
of star formation have minimal impact on the global 
galaxy color, the residual star formation had been un- 
constrained. These results will therefore provide useful 
input for models that attempt to measure the ages of 
the stellar populations of galaxies based on their spec- 
tral features. In particular, it would be inte resting to see 
wheth er frosting models such as those of iTrager et al.l 
(200.3), when using these SFHs, are able to reconcile the 
apparent discrepancy between the young ages inferred 
by the spectra compared to the old ages at which the 
majority of stars must have been formed. 

Future work includes a similar study of other nearby 
ellipticals and also lenticular galaxies, searching for in- 
dividual UV-bright young stars and open star clusters, 
to not only place tight constraints on their low-level star 
formation rates and histories, but to also address the role 
of environment and galaxy properties on star formation. 

This work has made use of BaSTI web tools. We 
thank Adriano Pietrinferni for providing stellar evolu- 
tionary tracks for the relevant WFC3 UVIS filters that 
were not publicly available. We also thank Jeremy Bailin, 
Eric Bell, and Sally Oey for providing valuable comments 
and assistance. Support for program # 11583 was pro- 
vided by NASA through a grant from the Space Telscope 
Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc,. under 
NASA contract NAS 5-26555. The National Radio As- 
tronomy Observatory is a facility of the National Sci- 
ence Foundation operated under cooperative agreement 
by Associated Universities, Inc. 
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